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Abstract 

We discuss the problem of finding "marginal" distributions within different tomographic 
approaches to quantum state measurement, and we establish analytical connections among 
them. 

PACS number(s): 03.65.-w, 42.50.-p 



1 Introduction 

Recently the problem of measuring quantum states in quantum optics and quantum mechanics 
found experimental realization |l|, 0, ^ Different procedures for measuring quantum states were 
suggested, i.e., to measure an observable and from the experimental result of the observations 
to reconstruct the density matrix of the quantum system in terms, for example, of its Wigner 
function or its quasidistribution function like the Husimi Q-function and the Glauber- 
Sudarshan P-distribution |^, ^. 
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The special interest to the problem of measuring the quantum states is related to the ex- 
perimental realization of specific nonclassical quantum states |10|- Such nonclasssical states, 



as one-mode even and odd coherent states |jTT| (or Schrodinger cat states) were created in high- 



quality cavity (see, for example, [|T2| ). It was also shown [|T^] that three-mode even and odd 
coherent states maight be created in such experiments. For an ion in a Paul trap, the even and 
odd coherent states may be realized 0, |14[ as well as the new type of nonclassical states like 



nonlinear coherent states 16, 17 



Thus, the problem of measuring such states, i.e., reconstructing their density matrix in any 
representation, is actual. For trapped ion several methods were indeed suggested for such a 
reconstruction [|1^]. Instead to determine the state of a radiation field there exists a dominant 
method, called optical tomography method which is based on a relation of the Wigner function 



to the distribution function of a homodyne observable [p!9i [20(1 , which may be derived from the 
general scheme of s-ordered quasidistributions |2ll . 



Recently, the optical tomography method was extended into symplectic tomography [^], 
in which the basic observable is the generic squeezed and rotated quadrature. It was used to 
find classical-like description of the quantum dynamic evolution If in optical tomography 

the Wigner function is reconstructed from measurable marginal distribution of the homodyne 
observable through Radon transform, in the symplectic tomography the Fourier transform is 
applied to reconstruct the density matrix by measuring the marginal distribution of the squeezed 
and rotated quadrature. In this context, the symplectic tomography is similar to the field-strength 
method of Ref. The symplectic tomography was generalized for multimode case as well [p5|]. 

Recently, it was also suggested the photon number tomography method |2^, ^7\, ^ , for which 
the measurable observable is the discrete photon number of the measurable field superposed with 
the local oscillator field by scanning its complex amplitude. The photon number distribution, de- 
pending on the complex amplitude of the local oscillator field, may be converted into the density 
matrix of the measurable field state. The symplectic tomography of nonclassical states (namely, 
squeezed states and Schrodinger cat states) of an ion in a Paul trap was discussed in Ref. 



2 



The reconstruction of the density matrix of the squeezed vacuum, by using reproducible measure- 



ments |^0| of the distribution for discrete photon numbers, was recently reported in Ref. IQ. 

Thus, in the theoretical and experimental framework of measuring quantum states there ex- 
ist different methods yielding the density matrix from the corresponding measurable probability 
distributions. On the other hand, till now the relations among all these distributions were not 
clarified in details, to our knowledge. The aim of this work is to establish the invertible transfor- 
mations of the symplectic quadrature distribution to the photon number distribution as well as 
to study the connection of the symplectic tomography scheme to the optical tomography. Some 
aspects of this problem were discussed in Ref. [^. The relation of the Radon transform to the 
symplectic tomography procedure was discussed in Ref. [^. By finding the explicit relations 
among the different measurable probability distributions one is aimed to make prediction of mea- 
surements for different observables, from measuruments, for example, of the marginal distribution 
of the generic squeezed and rotated quadrature in the symplectic tomography scheme. 

The paper is organized as follows. In Section 2, we review the properties of the marginal 
distributions extending the concept of " marginalizations" from the optical tomography to the 
symplectic, and photon number tomography. In Section 3, the connections among the various 
"marginal" distributions are given. The main results obtained are summarized in the Conclusion. 

2 Marginal Distributions 



By refering to the standard definitions given in the literature ^ , by "marginalization" one 
should intend a line integral in the phase space {q,p} of the Wigner function W{q,p), i.e. 

dqdp . 



W[X, 



-—-W{q,p)5{x-cos9q-smep), (1) 



where 9 is the angle orientation of the line. On the other hand, the marginal distribution w{x, 9) 
can be considered as well as the Fourier transform of the characteristic function xi^) 

w{x,9) = |^x(fc)e-^'^ (2) 
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with 

The marginal distribution represents a true measurable probability, at least in the field of optics, 
by means of homodyne detection methods. Furthermore, recently has been shown [|T9| that from 
a collection of marginal distributions it is possible to recover the phase space pseudo-distributions 
(tomography) or, even better, the density operator |Q representing the quantum state of the 
system under study. 

Our purpouse is now to generalize the previous definition of the marginalization in order to in- 
clude all measurable probabilities related to observables which allow a tomographic reconstruction 
of the quantum state . 

In fact, the tomographic principle may be generalized as follow: given a density operator p 
and a group element Q (or to better say an operator belonging to representation of a group, which 
acts in the space of quantum states), one can create different types of tomography if, by knowing 
the matrix elements 

w{y,g) = {y\gpg-'\y) (4) 

from measurements, is able to invert the formula expressing the density operator in terms of 
the above "marginal" distribution w {y,Q) , where y may denote either continuous or discrete 
eigenvalues. The positive distribution w {y,Q) is normalized 



Jw{y,g)dy = l. (5) 



It is to remark that, to simplify the terminology, from now on we refer indistincly to any positive, 
measurable, and normalized probability distribution like of Eq. (|^) as marginal distribution. 

To get the density operator from a marginal distribution an inversion procedure should be 
employed, and one can use the properties of summation or integration over group parameters Q. 
Thus, the marginal distribution of the observable y depends on extra parameters determining the 
operator Q. It is worth noting that the operator Q may belong not only to a group but also to 
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other algebraic construction, for example, to quantum group. The only problem is mathematical 
one, to make the inversion, and/or physical one, to realize the transformation Q the in laboratory. 

Along this line the "symplectic tomography" has been introduced in Ref. ||2^ by means of the 
transformation 

X = fiq + up (6) 

with the real parameters fi, v generalizing the homodyne rotation. In this case the marginal 
distribution becomes 

/dc^ dp 
^^-p p)6{x- fiq- up) , (7) 

while the density operator can be written in the invariant form as 

p = J dx dfidu w{x, fi,i')Kf^^^ , (8) 

with the kernel operator 

^ _]_^2^-izx^-z{u-itJ.)dyV2^z{u+itJ.)a/V2^~z^{n'^+iy^)/i ^g^ 

Here z represents the variable conjugate to x by the Fourier transform. The fact that K^,u 
depends on the z variable as well (i.e., each Fourier component gives a selfconsistent kernel) shows 
the overcompleteness of information achievable by measuring the observable of Eq. (^. Due to 
this, one could set 2; = 1 in (|). 

Beside that, a "photon number tomography" was developed in Ref. |^ parallely to Refs. |^ 
P7[ . It is based on the possibility to measure the number of photons in the state displaced over 
the phase space. Hence, the marginal distribution will be 

w(n, a) = TT{D{a)pD-\a)\n){n\} = {n\D{a)pD-\a)\n) , (10) 



and it also corresponds to a "propensity" obtained via filtering the original state described 



by p with the filter in a Fock state. The propensities have been also studied in connection with 
a quantum state reconstruction [^. In some sense the probability w{n,a) differs from the other 
distribution w(x, /i, z/) since it is intrinsically a phase space distribution (in the case of n = it 



corresponds to the Husimi-Q function). However, for convenience, we still refer to it as a marginal 
distribution in the sense previously explained. 



The invariant expression for the density operator in terms of the distribution (|l^) is 



P = y2 [ — w {n, a) ks{n,a) , [U] 



with the kernel given by 

KJn. a) = 

1-s Vs-1 
The marginal distribution is normalized as 



2 / s -I- 1 \ " 

K,in,a) = - (—-) Ti-a,-s). (12) 



^ w (n, a) = 1 . (13) 

n=0 

The operator T represents the complex Fourier transform of the s-ordered displacement operator 
D{^,s) = Z)(Oe'l«l'/^ which can also be written as IBll 



T{a,s) = ^—D{a)(^—^) D-\a) . (14) 



1-s ^ ' Vs-1 

3 Connection among marginal distributions 

Here we consider the possibility to express the above marginal distributions as function one of 
other. 

First, since the operator (^) is proportional to the displacement operator generating the coher- 
ent state from the vacuum, one can use the known matrix elements of the displacement operator 
in the Fock basis, expressed in terms of Laugerre polinomials to get the density matrix in 



the photon number basis as a convolution of the magrinal distribution w (a;, /i, v) 

[ZTn(n~m)/2 /7y2 i ,,2\ 

{m\p\n) = \ — / e*^-(^ )/^w(x, /i, i/)(i/ - 2/i)"^-"L"^-" ——^ dxdfxdu; m > n . 

\ ml 2n J \ 2 J 

(15) 

For m < n, we use the hermiticity property of the density matrix (m|p|n) = {n\p\m)*. The 
diagonal matrix elements of the density matrix in the Fock basis give the photon distribution in 
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terms of the marginal distribution 



{n\p\n) = P{n) = — / e*^-^'^ w{x, /i, z/) L„ dx dfi dv . (16) 

2tx J \ 2 j 

Let us now rewrite Eq. (0) as 

w{n, a) = Y,iji\D{a)\k){k\p\l)iJ\D~\a)\n) , (17) 

k,l 

then we can use the density matrix elements of Eq. (|T5|) and the number representation of the 



displacement operator given in Ref. [21|, to get 



+ z/^ aiv + iii) a*(i' — ijj) 
w [n, a) = — I w [x, /i, u) exp <tx h 



xL„ ( — j dxdpdv. (18) 

For a = , Eq. (|18D gives w {n, 0) = P{n) , where P{n) is the photon distribution (0). 

To derive the inverse relation, we start with an oscillator in the thermal equilibrium state. 
One can easily obtain the quadrature marginal distribution for the thermal equilibrium state of 
the harmonic oscillator, which is interacting with a heat bath at dimensionless temperature T, in 
the form 



-1/2 

Wt{x, h, u) = vr (/i^ + z/^) coth [{2T)^'^'^ 



exp 



x'^ 



(19) 



(/i2 + z/2) coth[(2T)-i] 

It means that for the non normalized state density operator of the harmonic oscillator of the form 

Prion = exp I — ^ I (20) 
we have the corresponding non normalized marginal distribution of the form (pnon ^non) 

If the density operator corresponds to the shifted equilibrium position in the driven oscillator 
phase space, i.e., 

Pshift = D{a) p^onD{-a) = exp ( j;^ ^ j , (22) 



we have the non normahzed marginal distribution of the driven oscillator 



Wa,s[X^ /i, U) = Wrp \X — II 



a + a* 



a — a 



V2 



— V 



72 



/i, V 



(23) 



One can see that the temperature T is 



T 



In 



s-\ 



s + 1 



(24) 



and the operator [(1 — s)/2] T(a, s) has exactly the form of the operator (|22D . It means that the 
non normalized marginal distribution corresponding to the operator T(q;, s) is 



w„,s(x, /i, v) 



1 - s 



TT + z/2) coth [(2T) 



-ii 



-1/2 



X exp < 



2 l-exp(-l/T) 

(x -[ii{a + a*) /v^] - [z/ (a - a*) /{i ^2)] ) 



, 2 ^ 



(/i2 + Z/2) coth [(2r)-l] 



(25) 



where the temperature T is given by (p^). Expression (pSj) gives the connection of the marginal 
distribution w (x, /i, i/) with the photon number marginal distribution w (n, a). In fact, by taking 
the expectation value of both side of Eq. (0) over the eigenket \x) of the quadrature operator 
(|§), one obtains 

(fa ^ 2 /s + l^'" 

m=0 

where 



w (x, /i, z/) = / > w (m, a) .^(x, /x, z/) 

J vr .ZTn i — s Vs — i/ 



(26) 



1 + S 



1 



TT (/i^ + z/2) coth [(2T) 



-ii 



-1/2 



X exp 



2 l-exp(-l/T) 

(x+ [/i(a + a*)/v/2] + [z/ (a - a*)/(z ^2)]) 
(/i2 + z/2) coth [(2r)-i] 



(27) 



and the "temperature" T is given by 



with s — > — s , i.e., 

s + 1^ 



T= In 



s - 1 



(2^ 



The equation (p6D , along with (pTl) and (^), is the inverse relation for (plsl) . Thus, measuring 
the quadrature marginal distribution w (x, /i, z/) in homodyne experiments one can predict the 



results of measurement of the photon marginal distribution and, vice versa, measuring the photon 
marginal distribution one can predict the results of measurements of the quadrature marginal 
distribution. 

Let us now address the question how to find the rotated quadrature marginal distribution 
w (x, 6) of the optical tomography scheme if one knows the squeezed and rotated quadrature 
distribution w (x, /i, v) of symplectic tomography. The answer is obvious, namely, 

w (x, cos 0, sin 9) = w (x, 9) . (29) 

Formula (p9|) may be rewritten in the integral form 

w {x, 9) = w (x, fi, u) 5 {fi — cos 9) 6 {u — sin 9) dfi du . (30) 



So, given the function w (x, fi, v) one gets the homodyne marginal distribution w (x, 9) . The 
inverse problem may also be solved, namely, how to find the marginal distribution of the symplectic 
tomography w (x, /i, v) , if one knows the homodyne marginal distribution w (x, 9) . It means that 
knowing the distribution function of two variables one has to reconstruct the distribution function 
of three variables. 

To make the reconstruction, we first decompose the periodic function w (x, 9) into the Fourier 
series 

oo 

w (x, 9) = [cn (x) COS n9 + dn (x) sin n9] . (31) 

n=0 

Then, in view of the known espressions for cos n9 , sin n9 as polynomials in sin 9 , cos 9 , we make 
in (|3l|) the replacement 

cos 9 ^ fi] sin 9 u . 
Thus, we get the expression of the marginal distribution w (x, /i, z/) in the form 



1 /■27r 1 °° r /■2i" 

/ w{x, 9) d9 + — y + iuf / e"*"^ w (x, 9) d9 + c.c. 

Jo 271 Jo 



W [X, /i, Z/J = — 

ZTT JO 



(32) 



Finally, by knowing the marginal distribution of the optical tomography scheme we directly obtain 
the marginal distribution of symplectic tomography. 
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Using (^) one can find the connection of the photon number tomography to the optical 
tomography formahsm. The marginal distribution of the rotated quadrature is expressed in terms 
of the photon number marginal distribution as follows 

(Pa ^ 2 + 



w 



where 



w^a,~s {x, COS 6, sin 



1 + s 1 



TT coth [(2r) 



-ii 



-1/2 



X exp 



2 1 - exp (-1/T) 

(x+ [cos ^(a + a*)/V2] + [sin 6 {a - a*)/{i ^2)]) 



. 2 ^ 



coth [(2T)-i] 



(34) 

and the "temperature" T is given by (p^). 

4 Conclusion 

We have provided analytical relations which connect the "marginal" distributions of different 
tomographic processes. Physically they will give predictions on how to make cross checking of the 
state measurements or how, by using one sort of measurement, to reconstruct another sort, i.e., 
by measuring an observable one gets information on the distribution of other observables. 

References 

[1] D. T. Smithey, M. Beck, M. G. Raymer, and A. Faridani, Phys. Rev. Lett. 70, 1244 (1993). 

[2] D. M. Meekhof, G. Monroe, B. E. King, W. M. Itano, and D. J. Wineland, Phys. Rev. Lett. 
76, 1796 (1996). 

[3] I. J. Dunn, I. A. Walmsley, and C. Mukamel, Phys. Rev. Lett. 74, 884 (1995). 

[4] S. Schiller, G. Breitenbach, S. F. Pereira, T. MuUer, and J. Mlynek, Phys. Rev. Lett. 77, 
2933 (1996). 

10 



[5] E. P. Wigner, Phys. Rev. 40, 749 (1932); E. P. Wigner, in Perspectives in Quantum Theory, 
W. Yourgrau and A. van der Merwe eds. (Dover, New York, 1979). 

[6] K. Husimi, Proc. Phys. Math. Soc. Jpn. 23, 264 (1940). 

[7] R. J. Glauber, Phys. Rev. Lett. 10, 87 (1963). 

[8] E. C. G. Sudarshan, Phys. Rev. Lett. 10, 277 (1963). 

[9] V. Buzek and P. L. Knight, in Progress in Optics, Ed. E. Wolf (Elsevier, Amsterdam, 1995), 
Vol. XXXIV, p. 1. 

[10] V. Buzek, C. Adam, and G. Drobny, Phys. Rev. A 54, 804 (1996); Ann. Phys. 245, 37 (1996). 
[11] V. V. Dodonov, I. A. Malkin, and V. L Man'ko, Physica 72, 597 (1974). 
[12] S. Haroche, Nuovo Cim. B 110, 545 (1995). 
[13] C. C. Gerry, Phys. Rev. A 54, R2529 (1996). 

[14] R. L. de Matos Filho and W. Vogel, Phys. Rev. Lett. 76, 608 (1996). 
[15] M. M. Nieto, Phys. Lett. A 219, 180 (1996). 

[16] V. I. Man'ko, G. Marmo, E. C. G. Sudarshan, and F. Zaccaria, in Proceedings of the Fourth 
Wigner Symposium (Guadalajara, Mexico, July 1995), eds. N.M. Atakishiyev, T.H. Seligman, 
and K.-B. Wolf (World Scientific, Singapore, 1996), p. 421; "f-Oscillators and nonlinear 
coherent states," Los Alamos Report No. |quant-ph/961200"^ ; Phys. Scripta (to appear). 



[17] R. L. de Matos Filho and W. Vogel, Phys. Rev. A 54, 4560 (1996). 

[18] see e.g.: S. Wallentowitz and W. Vogel, Phys. Rev. Lett. 75, 2932 (1995); J. F. Poyatos, 
R. Walser, J. J. Cirac and P. Zoller, Phys. Rev. A 53, R1966 (1996); C. D'Helon and G. J. 
Milburn, Phys. Rev. A 54, R25 (1996); P. J. Bardroff, C. Leichte, G. Schrade, and W. P. 
Schleich, Phys. Rev. Letts. 77, 2198 (1996). 

11 



[19] K. Vogel and H. Risken, Phys. Rev. A 40, 2847 (1989). 

[20] J. Bertrand and R Bertrand, Found. Phys. 17, 397 (1987). 

[21] K. E. Cahill and R. J. Glauber, Phys. Rev. 177, 1857; 1882 (1969). 

[22] S. Mancini, V. I. Man'ko, and P. Tombesi, Quantum Semiclass. Opt. 7, 615 (1995). 

[23] S. Mancini, V. I. Man'ko, and P. Tombesi, Phys. Lett. A 213, 1 (1996). 

[24] H. Kiihn, D. G. Welsch, and W. Vogel, J. Mod. Opt. 41, 1607 (1994). 

[25] G. M. D'Ariano, S. Mancini, V. I. Man'ko, and P. Tombesi, Quantum Semiclass. Opt. 8, 1017 
(1996). 

[26] S. Wallentowitz and W. Vogel, Phys. Rev. A 53, 4528 (1996). 

[27] K. Banaszek and K. Wodkiewicz, Phys. Rev. Lett. 76, 4344 (1996). 

[28] S. Mancini, V. I. Man'ko, and P. Tombesi, Europhys. Lett. 37, 79 (1997). 

[29] O. V. Man'ko, J. Russ. Laser Research (Plenum Press) 17, 439 (1996). 

[30] L. D. Landau and R. Peierls, Z. Phys. 69, 56 (1931). 

[31] S. Mancini, V. I. Man'ko, and P. Tombesi, "Different Realizations of Tomographic Principle 
in Quantum State Measurement", J. Mod. Opt. (submitted). 

[32] A. Wiincshe, "Radon Transform and Pattern Functions in Quantum Tomography", J. Mod. 
Opt. (to appear). 

[33] G. M. D'Ariano, U. Leonhardt, and H. Paul, Phys. Rev. A 52, R1801 (1995). 
[34] K. Wodkiewicz, Phys. Rev. Lett. 52, 1064 (1984). 

[35] V. Buzek, C. H. Keitel and P. L. Knight, Phys. Rev. A 51, 2575 (1995); ibid. 52, 2594 (1995). 

12 



